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Abstract These lecture notes on 2D growth processes are divided in two parts. 

The first part is a non-technical introduction to stochastic Loewner 
evolutions (SLEs). Their relationship with 2D critical interfaces is il- 
lustrated using numerical simulations. Schramm's argument mapping 
conformally invariant interfaces to SLEs is explained. The second part 
is a more detailed introduction to the mathematically challenging prob- 
lems of 2D growth processes such as Laplacian growth, diffusion limited 
aggregation (DLA), etc. Their description in terms of dynamical con- 
formal maps, with discrete or continuous time evolution, is recalled. We 
end with a conjecture based on possible dendritic anomalies which, if 
true, would imply that the Hele-Shaw problem and DLA are in different 
universality classes. 

Growth phenomena are ubiquitous in the physical world at many 
scales, from crystals to plants to dunes and larger. They can be studied 
in many frameworks, deterministic of probabilistic, in discrete or con- 
tinuous space and time. Understanding growth is usually a very difficult 
task. This is true even in two dimensions, the case we concentrate on in 
these notes. 

Yet two dimensions is a highly favorable situation because it allows 
to make use of the power of complex analysis in one variable. In many 
interesting cases, the growing object in two dimensions can be seen as 
a domain, i.e. a contractible open subset of the Riemann sphere (the 
complex plane with a point at infinity added). A deep theorem of Rie- 
mann asserts that such a domain, whatever complicated and fancy, is 
conformally equivalent to a simple reference domain, which is usually 
taken as the upper-half plane or the unit disk. This conformal equiva- 
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lence is unique once an appropriate normalization, which may depend 
on the growth problem at hand, has been chosen. Cauchy's theorem 
allows to write down an integral representation for the conformal map 
as an integral along the boundary of the reference domain, involving a 
density. This density is time dependent. Then the time derivative of 
the conformal map has an analogous representation and a nice way to 
specify the growth rule is often directly on this density. This leads to 
the concept of Loewner chains, which is the central theme of these notes. 
We shall illustrate Loewner chains in several situations. 

Our aim is to give a pedagogical introduction to a beautiful subject. 
We wanted to show that it leads to many basic mathematical structures 
whose appearance in the growth context is not so easy to foresee, like 
Brownian motion, integrablc systems and anomalies to mention just a 
few. We have also tried to stress that some growth processes have rules 
which are easy to simulate on the computer. A few minutes of CPU 
are enough to get an idea of the shape of the growing patterns, to be 
convinced that something interesting and non trivial is going on, and 
even sometimes to get an idea of fractal dimensions. This is of course 
not to be compared with serious large scale simulations, but it is a good 
illustration of the big contrast between simple rules, complex patterns 
and involved mathematical structures. However, other growth models, 
and among those some have been conjectured to be equivalent to simple 
ones, have resisted until recently to precise numerical calculations due 
to instabilities. 

To avoid any confusion, let us stress that being able to describe a 
growth process using tools from complex analysis and conformal geome- 
try does not mean that the growth process itself is conformally invariant 
at all. Conformal invariance of the growth process itself puts rather 
drastic conditions on the density that appears in the Loewner chain. 

This is illustrated by the first part of these notes, which deals with 
conformally invariant interfaces and their relation to stochastic Loewner 
evolutions. This part is an elaboration of the the main points developed 
during the lectures. 

The second part is an introduction to a larger class of processes de- 
scribing the growth of possibly random fractal planar domains and a 
review of some of their basic properties. Due to lack of time, this part 
was not presented during the lectures. 

The study of the continuum limit of non-intersecting curves on the 
lattice has been a subject of lasting interest both in mathematics and in 
physics. A famous example is given by self-avoiding random walks. The 
motivation comes from combinatorics, but also from statistical mechan- 
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ics. Two dimensions is most interesting because a non-intersecting curve 
is the boundary between two domains and can very often be interpreted 
as an interface separating two coexisting phases. 

At a critical point and for short range interactions, such interfaces are 
expected to be conformally invariant. The argument for that was given 
two decades ago in the seminal paper on conformal field theory [1]. The 
rough idea is the following. At a critical point, a system becomes scale 
invariant. If the interactions on the lattice are short range, the model 
is described in the continuum limit by a local field theory and scale 
invariancc implies that the stress tensor is tracclcss. In two dimensions 
this is enough to ensure that the theory transforms simply -no dynamics 
is involved, only pure kinematics- when the domain where it is defined 
is changed by a conformal transformation. 

The local fields arc classified by representations of the infinite dimen- 
sional Virasoro algebra and this dictates the way correlation functions 
transform. This has led to a tremendous accumulation of exact results 
using conformal field theory (CFT). A situation that is well under control 
is that of unitary minimal models. The Hilbert space of the system splits 
as a finite sum of representations of the Virasoro algebra, each associated 
to a (local) primary field, and the corresponding correlation functions 
can be described rather explicitly. The study of non-local objects like 
interfaces at criticality has not seen such a systematic development and 
only isolated though highly artful results [2, 3, 12] have been discovered 
using conformal field theory techniques. 

Non local objects like interfaces are not classified by representations 
of the Virasoro algebra but the reasoning that led O. Schramm to the 
crucial breakthrough [4], i.e. the definition of stochastic Loewner evolu- 
tions (SLEs), rests on a fairly obvious but cleverly exploited statement of 
what conformal invariance means for an interface. Surprisingly it allows 
to turn this problem into growth problem, something which looks natu- 
ral only a posteriori. For some years, probabilistic techniques have been 
applied to interfaces, leading to a systematic understanding that was 
lacking on the CFT side. A sample, surely biased by our ignorance, can 
be found in refs.[5, 6, 7, 8, 9]. There is now a satisfactory understand- 
ing of interfaces in the continuum limit. However, from a mathematical 
viewpoint, giving proofs that a discrete interface on the lattice has a 
conformally invariant limit remains a hard challenge and only a handful 
of cases has been settled up to now. 

There is now a good explanation of the -initially mysterious- rela- 
tionship between SLE and CFT. This was one of the main topics of the 
lectures. Though this aspect is mostly due to the work of the present 
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authors ([10], but see also [11]), we have decided to be very sketchy. The 
interested reader is referred to the hterature. 

Stochastic Loewner evolution is a simple but particularly interesting 
example of growth process for which the growth is local and continu- 
ous so that the resulting set is a curve without branching. Of course 
other examples have been studied in connection with 2d physical sys- 
tems. The motivations are sometimes very practical. For instance, is 
it efficient to put a pump in the center of oil film at the surface of the 
ocean to fight against pollution? The answer has to do with the Lapla- 
cian growth or Hele-Shaw problem [14]. The names diffusion limited 
aggregation [15], see Figure (1.1), and dielectric breakdown [16] speak 
for themselves. Various models have been invented, sometimes with less 
physical motivation, but in order to find more manageable growth pro- 
cesses. These include various models of iterated conformal maps [17], 
etc. See refs.[18] for recent reviews. As mentioned before, in most cases 
the shape of the growing domains is encoded in a uniformizing confor- 
mal map whose evolution describes the evolution of the domain. The 
dynamics can be either discrete or continuous in time, it can be either 
deterministic or stochastic. But the growth process is always described 
by a simple generalization of the Loewner equation called a Loewner 
chain. 




Figure 1.1 An example of DLA cluster obtained by iterating conformal maps. 

The general understanding of these Loewner chains is still embryonic, 
in strong contrast with the case of Loewner evolutions, which deal with 
local growth. For many problems the most basic question are unan- 
swered. For instance, as we shall explain below, Laplacian growth is at 
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the same time a completely integrable system and an ill-posed problem 
because it develops singularities in finite time. Hydrodynamics gives 
a natural regularization via the introduction of surface tension. It has 
been argued for some time that in the limit of vanishing surface tension 
one retrieves a model which is in the same universality class as diffu- 
sion limited aggregation. But the experimental and numerical evidence 
is inconclusive and there is no consensus. In fact, at the end of these 
notes we shall give an argument suggesting that they belong to different 
universality classes. This is very conjectural, but as should be clear from 
this introduction, non-local growth processes are certainly a source of 
interesting and challenging problems. 
Table of contents: 

- Section 1: Critical interfaces and stochastic processes. 

- Section 2: Loewner chains. 



The following pedagogical introduction to conformally invariant ran- 
dom curves breaks in two part. The first is a list of examples of critical 
interfaces on the lattice and some of their properties. The second is a 
presentation of O. Schramm's derivation of SLE. 

In this part the upper-half plane is used as a reference geometry. 



Let us start with three examples. Our aim is to explain their def- 
initions, to show a few samples and, in the first two cases, to give a 
numerical estimate of the fractal dimension. 

The first, loop-erased random walks, belongs to the realm of "pure 
combinatorics" . 

The second, percolation, is at the frontier between pure combinatorics 
and statistical mechanics because it arises very naturally as the domain 
boundary of a statistical mechanics model, but the Boltzmann weights 
are trivial. It is however a limiting case of a family of models with 
non trivial Boltzmann weights and, from that point of view, quanti- 
ties relevant for percolation arc obtained by exploring an infinitesimal 
neighborhood around the trivial Boltzmann weights. 

The third example, the Ising model, is deeply rooted in statistical 
mechanics. 



1. 



CRITICAL INTERFACES AND 
STOCHASTIC PROCESSES 



1.1 



THREE EXAMPLES AND A 
GENERALITY 
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We shall end this section by abstracting a crucial property of interfaces 
which allows to make an efficient use of conformal invariance and derive 
the Markov property for SLEs. 

1.1.1 Loop-erased random walks. This example is purely of 

combinatorial nature. A loop-erased random walk is a random walk with 
loops erased along as they appear. More formally, if Xq, Xi, - ■ ■ , is a 
finite sequence of abstract objects, we define the associated loop-erased 
sequence by the following recursive algorithm. 

Until all terms in the sequence are distinct. 

Step 1 Find the couple (/, m) with < I < m, such that the terms 
with indexes from to m — 1 are all distinct but the terms with indexes 
m and I coincide. 

Step 2 Remove the terms with indexes from Z + 1 to m, and shift the 
indexes larger than m by / — m to get a new sequence. 

Let us look at two examples. 
For the "month" sequence j, /, m, a, m, j, j, a, s, o, n, d, the first loop is 
m,a,m, whose removal leads to j, f,m,j,j,a,s,o,n,d, then j,f,m,j, 
leading to j, j, a, s, o, n, d, then j, j leading to j, a, s, o, n, d where all terms 
are distinct. 

For the "reverse month" sequence d,n,o,s,a,j,j,m,a,m,f,j, the first 
loop is leading to d,n,o,s,a,j,m,a,m,f,j, then a,j,m,a leading 
to d, n, o, s, a, m, f,j. 

This shows that the procedure is not "time-reversal" invariant. More- 
over, terms that are within a loop can survive: in the second example 
m, f, which stands in the j, m, a, m, f,j loop, survives because the first 
j is inside the loop a,j,m,a which is removed first. 

A loop-erased random walk is when this procedure is applied to a (two 
dimensional for our main interest) random walk. This is very easy to 
simulate. Fig. 1.2 represents a loop-erased walk of 200 steps. The thin 
lines build the shadow of the random walk (where shadow means that 
we do not keep track of the order and multiplicity of the visits) and the 
thick line is the corresponding loop-erased walk. The time asymmetry is 
clearly visible and allows to assert with little uncertainty that the walk 
starts on the left. 

To fit with the general SLE framework, let us restrict to loop-erased 
random walks in the upper-half plane. There are a few options for the 
choice of boundary conditions. 

A first choice is to consider reflecting boundary conditions on the real 
axis for the random walk. 

Another choice is annihilating boundary conditions: if the random walk 



Loewner Chains. 7 




Figure 1.2 A loop-erased random walk with its shadow. 

hits the real axis, one forgets everything and starts anew at the origin. 
Why this is a natural boundary condition has to wait until section 1.1.4. 

Due to the fact that on a two-dimensional lattice a random walk is 
recurrent (with probability one it visits any site infinitely many times), 
massive rearrangement occur with probability one. This means that if 
one looks at the loop-erased random walk associated to a given random 
walk, it does not have a limit in any sense when the size of the random 
walk goes to infinity. Let us illustrate this point. The samples in fig.1.3 
were obtained with reflecting boundary conditions. It takes 12697 ran- 
dom walk steps to build a loop-erased walk of length 633, but step 12698 
of the random walk closes a long loop, and then the first occurrence of 
a loop-erased walk of length 634 is after 34066 random walk steps. Ob- 
serve that in the mean time most of the initial steps of the loop-erased 
walk have been reorganized. 

However, simulations are possible because when the length of the 
random walk tends to infinity, so does the maximal length of the corre- 
sponding loop-erased walk with probability one. 

Though annihilating boundary conditions lead to remove even more 

parts of the random walk than the reflecting ones, the corresponding 
process can be arranged (conditioned in probabilistic jargon) to solve 
the problem of convergence as follows. 

Instead of stopping the process when the loop-erased walk has reached 
a given length, one can stop it when it reaches a certain altitude, say n, 
along the y-axis. Whatever the corresponding random walk has been, 
the only thing that matters is the last part of it, connecting the ori- 
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Figure 1.3 On the left: a large loop is about to be created. On the right: the massive 
rearrangement to go from 633 to 634 steps. 

gin to altitude n without returning to altitude 0. Moreover, the first 
time the loop-erased walk reaches altitude n is exactly the first time the 
random walk reaches altitude n. Now a small miracle happens: if a Id 
symmetric random walk is conditioned to reach altitude n before it hits 
the origin again, the resulting walk still has the Markov property. It 
is a discrete equivalent to the 3d Bessel process (a Bcsscl process de- 
scribes the norm of a Brownian motion, however no knowledge of Bessel 
processes is needed here, we just borrow the name). When at site m, 
< m < n, the probability to go to m it 1 is (1 ± l/m)/2, independently 
of all previous steps. Observe that there is no n dependence so that 
we can forget about n, i.e. let it go to infinity. The discrete 3d Bessel 
process is not recurrent and tends to infinity with probability one: for 
any altitude I there is with probability one a time after which the dis- 
crete 3d Bcsscl process remains above / for ever. Henceforth, we choose 
to simulate a symmetric simple random walk along the x axis and the 
discrete 3d Bessel process along the y-axis and we erase the loops of this 
new process. This leads to the convergence of the loop-erased walk and 
numerically to a more economical simulation. 

Fig. 1.4 is a simulation of about 10"" steps, both for reflecting and anni- 
hilating boundary conditions. At first glance, one observes in both cases 
similar simple (no multiple points) but irregular curves with possibly 
fractal behavior. 

To estimate the Hausdorff dimensions in both cases, we have gener- 
ated samples of random walks, erased the loops and made the statistics 
of the number of steps S of the resulting walks compared to a typical 
length L (end-to-end distance for reflecting boundary conditions, max- 
imal altitude for annihilating boundary conditions). In both cases, one 
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Figure I.4 A sample of the loop-erased random walk for the two boundary conditions. 

observes that S (X and a modest numerical effort (a few hours of 
CPU) leads to 5 = 1.25 ± .01. This is an indication that the boundary 
conditions do not change the universality class. 

To get an idea of how small the finite size corrections are, observe 
fig.1.5. The altitude was sampled from 2'^ to 2^^. The best fit gives a 
slope 1.2496 and the first two points already give 1.2403. 




2 4 6 8 10 



Figure 1.5 The logarithm of the average length of the loop-erased random walk 
versus logarithm of the maximum altitude. The numerical results are the circles, the 
straight line is the linear regression, the error bars are shown. 
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As recalled in the introduction, it is believed on the basis of intuitive 
arguments that in two dimensions scale invariance implies conformal in- 
variance, providing there are no long range interactions. What does this 
absence of long range interactions mean for loop-erased random walks? 
Clearly along the loop-erased walk there are long range correlations, if 
only because a loop-erased random walk cannot cross itself. However, 
the relevant feature for the intuitive argument is that, in the underlying 
2d physical space, interactions are indeed short range. At each time 
step, the increment of the underlying random walk is independent of the 
rest of the walk, and the formation of a loop to be removed is known 
from data at the present position of the random walk. 

iFvom the analytical viewpoint, the loop-erased random walk is one 
of the few systems that has been proved to have a conformally invariant 
distribution in the continuum limit, the fractal dimension being exactly 
5/4. A naive idea to get directly a continuum limit representation of 
loop-erased walks would be to remove the loops from a Brownian motion. 
This turns out to be impossible due to the proliferation of overlapping 
loops of small scale. However, the SLE2 process, to be defined later, 
gives a direct definition. In fact, it is the consideration of loop-erased 
random walks that led Schramm [4] to propose SLE as a description of 
interfaces. 

1.1.2 Percolation. To define a random interface for percolation, 
one possibility is to pave the upper-half plane with regular hexagons. 
Each one is, independently of the others, occupied or empty with prob- 
ability 1/2, with the exception of the ones along the real axis, which 
are all empty on the positive real axis and occupied on the negative real 
axis, see fig.1.6. Then a continuous interface, starting at the origin and 
separating occupied and unoccupied hexagons, is uniquely defined. As 
before, this defines a simple curve on the lattice. 

The percolation interface has an obvious but very remarkable property 
that singles it out: locality. Locality means that the percolation interface 
does not depend on the distribution of occupied and empty sites away 
from itself. Equivalently, if D is a domain (i.e. in this discrete setting 
a connected and simply connected family of hexagons in the upper half 
plane) containing the origin, the law of the percolation interface in D 
before it hits the boundary of D for the first time is independent of the 
distribution of occupied/empty sites outside D. It is the same law as 
for what one would define as the percolation interface in D without ever 
mentioning the world outside D. 
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Figure 1.6 The definition of the percolation interface. 



This observation makes the percolation interface easy to simulate. 
Indeed, the construction of a percolation interface proceeds inductively, 
as shown on the fig. 1.7. 
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Figure 1.7 The percolation interface as a growth process. 



Fig. 1.8 shows a few samples of increasing size. 

The contrast with the previous example is rather striking. Even for 
small samples, the percolation interface makes many twists and turns. 
With the resolution of the figure, the percolation interface for large sam- 
ples does not look like a simple curve at all! The intuitive explanation 
why this does not occur for loop-erased walks is that if it comes close to 



12 




Figure 1.8 Samples of the percolation interface for increasing sizes. 



itself, then with large probability a few more steps of the random walk 
will close a loop. 

Getting an estimate of the fractal dimension of the percolation inter- 
face proceeds along lines similar to the first example. If the number of 
steps S of the percolation interface is compared to a typical length L 
one observes that S oc and a modest numerical effort (a few hours of 
CPU) leads to 5 = 1.75 ± .01. 

The percolation interface is also build by applying rules involving 
only a few nearby sites, and again it is expected on general non rigorous 
arguments that its scale invariance should imply its conformal invariance 
in the continuum limit. 

The percolation process is another one among the few systems that 
has been proved to have a conformally invariant distribution in the con- 
tinuum limit, the fractal dimension being exactly 7/4. As suggested 
by numerical simulations, the continuum limit does not describe simple 
curves but curves with a dense set of double points, and in fact the -to 
be defined later- SLEg process describes not only the percolation inter- 
face but also the percolation hull [2], which is the complement of the set 
of points that can be joined to infinity by a continuous path that does 
not intersect the percolation interface. 
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1.1.3 The Ising model. Our next example makes also use of the 
same pictorial representation: it is the Ising model on the triangular 

lattice in the low temperature expansion. The spins arc fixed to be up 
on the left and down on the right of the origin along the real axis. The 
energy of a configuration is proportional to the length of the curves sep- 
arating up and down islands. The proportionality constant has to be 
adjusted carefully to lead to a critical system with long range correla- 
tions. This time, making accurate simulations is much more demanding. 
On the square lattice, the definition of the interface suffers from ambi- 
guities, but these become less relevant for larger sample sizes. 




Figure 1,9 A sample for the critical Ising model. The bottom line, where the spins 
are frozen -black on the right, white on the left- is not represented. Courtesy of J. 
Houdayer. 

Although there is no question that the fractal dimension of the Ising 

interface with the above boundary conditions is 11/8 and is described 
by -to be defined later- SLE3, a mathematical proof that a continuum 
limit distribution for the interface exists and is conformally invariant is 
still out of reach. 

We conclude this list of examples with two remarks. 

First, the pictorial representation by a gas of non-intersecting curves 
on the hexagonal lattice that we used in the second and the third example 
applies to another family of models, the 0{n) models, to which similar 
interface considerations would apply. 

Second, non branching interfaces, described on the lattice by simple 
curves, are not the generic situation. For instance the (5 = 3 states 
Potts model, three different phases coexist and the physical interfaces 
have branch points. 
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1.1.4 A generality. Up to now, we have discussed interfaces in 
-lattice approximations of- the upper half plane. Let us note that they 
make sense for more general domains. 

Start with the case of the percolation or of the Ising model. In the 
plane, take any connected and simply connected collection of hexagons, 
split it into an interior and a connected boundary, and split the boundary 
in two connected pieces in such a way that exactly two pairs of adjacent 
boundary hexagons, marked a and b, carry different colors, see fig. 1.10. 
Observe that we allow a cut which would correspond to the beginning 
of the interface. Hexagons separated by the cut are not to be counted 
as neighbors. 




Figure 1.10 A domain on which percolation and the Ising model can be defined. 

Then any coloring of the hexagons in the interior will lead to a well- 
defined interface. It will also lead to a well defined energy, and then the 
question of the distribution of the interface in such a domain is mean- 
ingful. For percolation the energy is independent of the configuration. 
For the Ising model, it is proportional to the length of the curves sep- 
arating up and down sites. Observe that -as the possible cut separates 
fixed spins- counting interactions across the cut or not in the energy 
just adds a constant to it, so that it has no influence on probabilities. 

For the loop-erased random case, the idea is similar. One takes a 
simply connected piece of the square lattice with two points a and b on 
the boundary, again allowing some cuts, see fig. 1.11. Consider all walks 
from a to 6 that do not touch the boundary except at a before the first 
step and at b after the last step and give each such walk of length I a 
weight 4~'. Observe that this choice is exactly the annihilating boundary 
condition: in the half plane geometry, a was the origin and b the point at 
infinity and, due to the infinite extension of the boundary, we had to go 
through a limiting procedure. Then erase the loops to get a probability 
distribution for loop-erased random walks from a to 6 in the domain. The 
probability for the simple symmetric random walk to hit the boundary 
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Figure 1.11 A domain on which loop-erased random walks can be defined. 

for the first time at h starting from a can be interpreted as the partition 
function for loop-erased walks. 

We can now go to the point we want to make, valid for all the above 
examples. For any of these, we use ^^(0,0,6) to denote the probability 
distribution for the interface ^[ab] from a to 6 in D. 

Suppose that we fix the beginning 7^^^] of a possible interface in do- 
main D, up to a certain point c. Then 1) we can consider the conditional 
distribution for the rest of the interface and 2) we can remove the be- 
ginning of the interface from the domain to create a new domain and 
consider the distribution of the interface in this new domain. This is 
illustrated on fig. 1.12. 




Figure 1.12 An illustration of situations 1) and 2) for the case of loop-erased walks. 
What is the distribution of the dotted curve in both situations ? 

We claim that the distributions defined in 1) and 2) coincide. For 
reasons to be explained in a moment, we call this property "locality at 
the interface" . In equations 



-P(D,a,b)( • \l[ac\) = -P(D\7[,,[,c,6)( • )• 
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It is obvious that these two probabilities are supported on the same 
set, namely simple curves along the edges of the lattice, going from c to 
6 in D \ 7[(ic[- Let us however note that for loop-erased random walks, 
annihilating boundary conditions are crucial. Reflecting boundary con- 
ditions clearly do not work, if only because the supports do not coincide 
in that case. 

For the case of loop-erased random walks, the argument for the equiv- 
alence of 1) and 2) goes as follows. Take any random walk (possibly 
with loops) Wo = a, Wi, ■ ■ ■ ,Wi = b that contributes to an interface j^^b] 
which is 7[(jc] followed by some 7[c6] • Let m be the largest index for which 
the walk visits c. Because the interface has to start with ^^^c] > the walk 
Wm = c,- ■ ■ ,Wi = b cannot cross ^y^acl again, so it is in fact a walk in 
D \ 7[ac[ from c to 6 leading to the interface j^cb] ■ The weight for the 
walk Wq = a, Wi, ■■■ ,Wi = b is A~\ i.e. simply the product of weights 
for the walks Wq = a,Wi, - ■ ■ , Wm = c and Wm = c = a, ■ ■ ■ ,Wi = b. 
Then a simple manipulation of weights leads directly to the announced 
result. 

For the case of percolation and the Ising model, in fact more is true: we 
can view -P(D,a,6) ^^t only as a probability distribution for the interface, 
but as the full probability distribution for the colors of the hexagons and 
still check the identity of 1) and 2). Again, the supports are the same for 
1) and 2), namely any configuration of the colors, except that the colors 
on both sides of j^ac] are fixed. For the case of percolation, the colors 
are independent of each other so the identity of 1) and 2) is clear. For 
the Ising model, the difference is that the conditional probabilities in 1) 
take into account the interactions between the colors along the interface, 
whereas the probability in 2) does not take into account the interactions 
between the spins along the cut left by the removal of the interface. 
However, as already mentioned above, the corresponding colors are fixed 
anyway, so the Boltzmann weights for the configurations that are in the 
support of 1) or 2) differ by a multiplicative constant, which disappears 
when probabilities are computed. 

This argument extends immediately to systems with only nearest 
neighbor interactions. They can be defined on any graph. If any sub- 
set of edges is chosen and the configuration at both end of each edge is 
frozen, it makes no difference for probabilities to consider the model on 
a new graph in which the frozen edges have been deleted. 

Instead of looking for further generalizations, we argue more heuristi- 
cally that the continuum limit for a system with short range interactions 
should satisfy locality at the interface. The use of this locality property 
-which, as should be amply evident, has nothing to do with conformal 
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invariance- together with the conformal invariance assumption is at the 
heart of O. Schramm's derivation of stochastic Loewner evolutions. 

1.2 O. SCHRAMM'S ARGUMENT 

We break the argument in several pieces. The heart of the prob- 
abilistic derivation establishes two properties of conformally invariant 
interfaces: 

■ the Markov property, 

■ the stationarity of increments, 

see sections 1.2.3, 1.2.2 and 1.2.5. The crucial results needed for the 
probabilistic part, in particular the basics of Loewner evolutions, are 
collected in 1.2.1 and 1.2.4. The last two subsections are not directly 
related to the argument. They collect some basic facts on the fractal 
properties of SLE interfaces and on the connection with CFT. 

1.2.1 Riemann's theorem and hulls. In this section we collect 
a few indispensable results used in the rest of the course. 

A domain is a non empty connected and simply connected open set 
strictly included in the complex plane C. Simple connectedness is a 
notion of purely topological nature which in two dimensions asserts es- 
sentially that a domain has no holes and is contractible. But it is a deep 
theorem of Riemann that two domains are always conformally equiv- 
alent, i.e. there is an invertible holomorphic map between them. For 
instance, the upper-half plane H is a domain. It is well known that it has 
a three dimensional Lie group of conformal automorphisms, PSL2(R,), 
that also acts on the boundary of H. There is a unique automorphism, 
possible followed by a transposition, that maps any triple of boundary 
points to any other triple of boundary points. By Riemann's theorem, 
this is also true for any other domain, at least if the boundary is not too 
wild. 

Riemann's theorem is used repeatedly in the rest of this course and 
is the starting point of many approaches to growth phenomena in two 
dimensions. 

For later use, we note that one can be a bit more explicit when the 
domain D differs only locally from the upper half plane H, that is if 
K = H \ D is bounded. Such a set K is called a hull. The real points in 
the closure of K in C form a compact set which we call Kr. Let / : H 
D be a conformal bijection. As the boundary of H is smooth, / has a 
continuous extension to R = R U oo, and /^^(R \ Kr,) is a non-empty 
open set in R with compact complement. We call the complement the 
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cut of /. By the Schwarz symmetry principle, defining f{z) = f{z) for 
SJm 2; < gives an analytic extension of / to the whole Riemann sphere 
minus the cut. Across the cut, / has a purely imaginary nonnegative 
discontinuity which we write as a Radon-Nikodym derivative djif /dx. 

One can use the PSL2 (R) automorphism group of H to ensure that 
/ is holomorphic at 00 and f{w) — w = 0{l/w) there. This is called the 
hydrodynamic normalization. It involves three conditions, so there is no 
further freedom left. We shall denote this special representative by /k, 
which is uniquely determined by K: any property of /k is an intrinsic 
property of K. 

Cauchy's theorem yields 

(1.1, 

A quantity that plays an important role in the sequel is 

a positive (unless K = 0) number called the capacity of K, which is 
such that /k(u^) = w — Ck/w + 0{1/w'^) at infinity. The usefulness of 
capacity stems from its good behavior under compositions: if K and K' 
are two hulls, K U /k(K') is a hull and 

Cku/k(K') = C'k + C'k', (1-2) 

as seen by straightforward expansion at infinity of /k o /k' ) the map as- 
sociated to KU/k(K'). In particular capacity is a continuous increasing 
function on hulls. 

Anticipating a little bit, let us note immediately that giving a dy- 
namical rule for the evolution of the finite positive measure d^Xf^ix) is 
a good way to define growth processes. 

1.2.2 Conformally invariant interfaces. Consider a domain 
D, with two distinct points on its boundary, which we call a and h. 

A simple curve, denoted by 7[ab]i from a to 6 in D is the image of a 
continuous one-to-one map 7 from the interval [0, -|-oo] to D U {a, 6} 
such that 7(0) = a, 7(00) = h and 7]afe[ = 7(]0, oo[) C D. Alternatively, 
a simple curve from a to 6 is an equivalence class of such maps under 
increasing reparametrizations. A point on it has no preferred coordinate 
but is has a past and a future. If c G D is an interior point, we use a 
similar definition for a simple curve ^^^c] from a to c in D. 
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Note that, apart from the fact that on the lattice we could use lattice 
length as a parameter along the interface, we have just rephrased in 
the continuum -but with the same notations- what we did before in a 
discrete setting. 

Our aim is to study conformally invariant probability measures on the 
set of simple curves from a to 6 in D. There is a purely kinematical step, 
which demands that if h is any conformal map that sends D to another 
domain /i(D), the measure for (/i(D), h{a), h{b)) should be the image by 
h of the measure for (D, a, b): 

P{D,a,b){l[ab] CU) = -P(/i(D),h(a),h(6))(^(7[afe]) C h{U)), 

where -P(D,a,6) (7[af<] denotes the probability for the curve j^^b] to 

remain in a subset U of D. See fig.1.13. 



This condition is natural and it is the one that conformal field theory 
suggests immediately. Let us note however that a totally different def- 
inition of conformal invariance is understood in the familiar statement 
"two dimensional Brownian motion is conformally invariant" . 

Observe that we could take any measure for (D, a, b) -well, with the 
invariance under the one parameter group of automorphisms that fixes 
(D, a, b)- and declare that the measure in h(D) is obtained by definition 
by the rule above. To make progress, we need to combine conformal 
invariance with locality at the interface. 

1.2.3 Markov property and stationarity of increments. This 
short section establishes the most crucial properties of conformally in- 
variant interfaces. 

Take c G D and let 7[ac] be a simple curve from a to c in D. Observe 
that D \ 7]ac] is a domain. To answer the question "if the beginning of 
the interface is fixed to be 7[ac]i what is the distribution of the rest 7|^jj 




a 



Figure 1.13 Conformal invariance for change of domain. 
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of the interface ?" we apply locahty at the interface to argue that this is 
exactly the distribution of the interface in D\7]ac]- We map this domain 
conformally to D via a map h-y sending 6 to 6 and c to a, see fig. 1.14. 




Figure I.I4 Conformal invariance for conditional probabilities. 



Conformal invariance implies that the image measure is the original 
measure in P(p^a,b)j where hj{'y'^^^) is an interface that has forgotten 
7[ac]- To summarize : 

^7('^[c&]) independent of j^g^^j (the Markov property) and has 
the same distribution as the original interface itself (stationarity 
of increments). 

1.2.4 Local growth and Loewner evolutions. The Markov 
and stationarity of increments property make it plain that to understand 
the distribution of the full interface, it is enough to understand the 
distribution of a small, or even infinitesimal, initial segment, and then 
glue segments via conformal maps. 

This calls for a description by differential equations. 

For this purpose, it is very convenient -even if by no means mandatory- 
to have a natural parametrization of interfaces. Using conformal in- 
variance, we can restrict ourselves to the situation when (D, a, b) = 
(H,0,oo). If 7[ooo] is a simple curve from to 00 in H, and c a point 
on it, we know that H \ tjqc] is a domain, and 7]oc] itself is a hull. 
We Tisc capacity as a parametrization and define a time parameter by 
2t{c) = C-^jo^j • The factor of 2 is just historical. The map i is a continu- 
ous increasing function of c varying from to 00 in R, when c moves from 
to 00 along 7[ooo]- The inverse map t ^ c{t) is well-defined and gives 
a parametrization of the curve. So we introduce an increasing family of 
hulls Kt,t G [0, oo[, by = 7]oc(t)] which has capacity 2t. Let ft = /kj 
be the conformal homeomorphism from H to H \ normalized to sat- 
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isfy ftiw) = w - 2t/w + 0(l/u)^) at infinity. Define : H \ Kj i-^ H to 
be the inverse of ft- Then gt{z) = z + 2t / z + 0{1 / z^) at infinity. 

To study the evolution of the family of hulls K^, fix e > and consider 
the hull Kg^t = gt{Vit+e \ Kf), which has capacity 2e by eq.(1.2). Define 
fe,t = /Ke,t • Then gt = fe,t o gt+e on H \ Kt+s- Using the representation 
of fKe,t terms of its discontinuity eq.(l.l), we obtain 



.X] 



gt+s — fe,t ° gt+e — gt+e — gt — TT' I '' — 

gt+e-x 

We introduce now the notion of local growth which is crucial for in- 
terfaces. When e is small, Kg^^ is a tiny piece of curve and the STipport 
of ci/x/jt is small and becomes a point when e goes to 0. Measures sup- 
ported at a point are S functions, so there is a point such that, as a 
measure, d^f^^/dx ~ 2e5{x — ^t) as e ^ 0^. If is a more general 
increasing family of hulls of capacity 2t, we say that the condition of 
local growth is satisfied if the above small e behavior holds. At first 
sight, it might seem that local growth is only true for curves, but this is 
not true. We shall give an example below. 

Letting e ^ 0"*", from the local growth condition, we infer the existence 
of a real function such that 

d9t^ . 2 

It is useful to look at this equation from a slightly different point of 
view, taking the function as the primary data. The solutions of this 
equation for a given function with initial condition go (z) = z is called 
a Loewner evolution. The image of by g^^^ is the tip of the curve at 
time t. 

The cases when the local growth condition is not satisfied are called 
Loewner chains, see below. Had we used another parametrization of the 
curve, the 2 in the numerator would be replaced by a positive function 
of the parameter along the curve. 

Informally, if is a growing curve, we expect that gt-^s{z) — gt.{z) 
describes an infinitesimal cut. This is confirmed by the explicit solution 
of eq.(1.3) for the trivial case = 0, which yields gt{z)^ = z^ + 4t, the 
branch to be chosen being such that at large z, gt{z) ~ z. This describes 
a growing segment along the imaginary axis. So intuitively, the simple 
pole in eq.(1.3) accounts for the existence of a cut and different functions 

account for the different shapes of curves. 

One can also solve the case when Kj is an arc of circle starting from 
the origin. It can be obtained from the trivial solution above by applying 
appropriate time dependent PSL2 (R) transformations both to gt and z 
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and then by a time change to recover the capacity parametrization. This 
is an illuminating exercise that we leave to the reader. Take an arc going 
from to 2R along a circle of radius R. It turns out that when the arc 
approaches the real axis to close a half disk, the function has a square 
root singularity oc — 2t. The capacity remains finite and goes 
to and the map itself has a limit (7^2 /2i^) = + ^ /{^ + which 
has swallowed the half disk without violating the local growth condition. 
One can start the growth process again. Making strings of such maps 
with various values of the radii is a simple way to construct growing 
families of hulls that are not curves and that nevertheless grow locally. 
Note that a square root singularity for is the marginal behavior: if 
is Holder of exponent > 1/2, Loewner evolution yields a simple curve. 

1.2.5 Stochastic (or Schramm) Loewner evolutions. If we 

sample locally growing hulls with a certain distribution, we get an as- 
sociated random process In the case of a conformally invariant dis- 
tribution, we have established two crucial properties: Markov property 
and stationarity of increments. To finish Schramm's argument leading 
to SLE, what remains is to see the implications of these properties on 
the distribution of 

The argument and expressions for the meaning of Markov property 
and stationarity of increments involved a map h that mapped the tip of 
the piece of interface to the initial marked point a and the final marked 
point h to itself. The map ht{z) = gt{z) — has the required property 
when the domain is the upper-half plane with and 00 as marked points. 
It behaves like ht{z) = z — + 2t/z + 0(l/z^) at infinity. We infer that 
for s > t, ht(K.s \ Kt) is independent of Kj/, t' < t (Markov property) 
and is distributed like a hull of capacity s — t = C/j^(Ks\Kt) (stationarity 
of increments). 

The hull determines the corresponding map h, so this can be rephrased 

as: hsoh~[^ (which uniformizes /it(Ks\Kj)) is independent of /if/, t' < t, 
and distributed like an As/igo/i^^ = z—{^s—£,t)+{s—t)/ z+0{l/ z"^) 
at infinity, the driving parameter for the process hg o h^^ is — Ct- To 
summarize: 



the Markov property and stationarity of increments for the inter- 
face lead to the familiar statement for the process ^j: for s > t, 
— is independent of ^t' , t' < t, (Markov property) and dis- 
tributed like a ^g-t (stationarity of increments). 

To conclude, a last physical input is needed: one demands that the 
interface does not branch, which means that at two nearby times the 
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growth is at nearby points. This imphes that is a continuous process, 
in the sense that it has continuous trajectories. 

One is now in position to apply a mathematical theorem: a Id Markov 
process with continuous trajectories and stationary increments is propor- 
tional to a Brownian motion. We conclude that there is a real positive 
number k such that = \/K,Bt for some normalized Brownian motion Bt 
with covariance E[i?si?j] = min(s,t). The same argument without im- 
posing that the time parametrization is given by the capacity of the hull 
would lead to the conclusion that the driving parameter is a continuous 
martingale, which is nothing but a Brownian motion after a possibly 
random time change. 

A solution of 

^(^) = 1 (14) 

dt ^ ' gt{z) - ^Bt ^ ■ ' 

is called a chordal Schramm-Loewner evolution of parameter k, in short 
a chordal SLE^^. The connection of this equation with interfaces relies 
mainly on conformal invariancc. But local growth, absence of branches 
and to a lower level locality at the interface also play a crucial role. 

1.2.6 Miscellanea on SLE. Up to now wc have only discussed 
the situation when the interface goes from to oo in the upper half 
plane, or more generally from one point on the boundary of a domain to 
another one. This is called chordal SLE. The argument can be extended 
easily to two other situations. The first is for an interface starting at 
a fixed point on the boundary of a domain and ending at a fixed point 
in the bulk. This is called radial SLE. For the second situation, three 
points are chosen on the boundary and the interface starts from the first 
point and ends at a random point between the two others. This is called 
dipolar SLE. Just as the upper-half plane is a convenient geometry for 
chordal SLE, a semi infinite cylinder is nice for radial SLE and an infinite 
strip is nice for dipolar SLE. 

For radial SLE on a cylinder of circumference tt, the equation reads 

^(z) = - fl5) 

dt ^ ' tan(c/t(z) - ^f^Bt) ^ ' ' 

and for dipolar SLE on a strip of width 7r/2, one has 

dgt , . ^ 2 

dt tanh((7t(z)- V^^O ^ ' 

The generalization of SLE to Riemann surfaces with moduli is easy for 
the annulus, but raises non-trivial problems in general, currently subject 
of active research. 
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The set of exact results obtained for SLE forms an impressive body of 
knowledge. We shall not mention the ones dealing with the explicit com- 
putation of certain crossing probabilities but wc list just a few "pictorial" 
properties with some comments. They (the properties and the com- 
ments) should be understood with the standard proviso "almost surely" 
or "with probability 1". 

We start with a surprisingly difficult result [4, 5, 6, 7]. 

■ Whatever the value of k, the pre-image of the driving parameter 
lim^_,^5^ 5(j~^(u)) is a continuous curve jt, called the SLE trace. 
The trace never crosses itself. This property is crucial if the trace 
is to be interpreted as a curve separating two phases. 

■ For K G [0, 4] the SLE trace is a simple curve. For k g]4, 8[, it has 
double points. For k G [8, oo[, it is space filling. 

■ The fractal dimension of the trace is 1 -|- k/S for k < 8 and 2 
for K>8. 

Using the formula for the dimension of the trace and confronting with the 
numerical simulations, it is plausible (actually, these are among the few 
cases for which a mathematical proof exists) that loop-erased random 
walks correspond to k = 2,c? = 5/4 and percolation to k = 6, d = 7/4. 
This is also compatible with the general shape of the numerical samples, 
which indicate that loop-erased random walks indeed lead to simple 
curves and that percolation doesn't. 

The hull K.t is by definition H\5^^(H). It has the following properties 

■ The hull is the complement of the connected component of oo 
in H\7]o,t]. 

■ For K £ [0,4], the SLE hull is a simple curve coinciding with the 
trace. For n g]4, oo[, the SLE hull has an nonempty connected and 
relatively dense interior. 

This may seem surprising at first sight. It is the sign that for k > 4, the 
drift i/k-Bj goes fast enough for the swallowing procedure to take place, 
as described in the closing arc example, but on all scales. 
This is summarized by fig. 1.15. 

For the many other properties known about SLE, the reader is invited 
to read the vast literature [4, 5, 6, 7, 8, 9]. 

1.3 SHORT REMARKS ON SLE AND CFT 

As already mentioned, it is by no means an easy task to prove that a 
discrete random curve converges to an SLE. For geometric models like 
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0<K!s4 4<K<8 K»8 



Figure 1.15 The phases of SLE. 

percolation, loop-erased random walks, self avoiding walks, it is often 
easier to compute the appropriate k using known special properties of 
the discrete model that are expected to survive in the continuum limit. 

For instance, one can formulate locality in the continuum. Showing 
that SLEk has the locality property only for k = 6 is a standard com- 
putation in stochastic Ito calculus, see e.g. [9]. 

For critical statistical mechanics models with non trivial Boltzmann 
weights, things are more complicated. The system contains local dy- 
namical degrees of freedom independently from the interface and, in the 
continuum limit, these local degrees of freedom are expected to be de- 
scribed by a conformal field theory (CFT). The most basic parameter of 
a CFT is its central charge c. Computing it (even non rigorously) from 
the lattice model can be a challenge comparable to the one of computing 

K. 

On the other hand, the relation between c and k can be worked out 
in general. The idea is elementary. From a technical viewpoint, it rests 
on assumptions similar to the ones needed in O. Schramm's argument. 

In the discrete setting, take O to be any observable. If one computes 
(0)d\7]„c] ill D \ 7]ac]) i-e. with part of the interface fixed, and then 
averages over 7]ac] one retrieves (0)d: 

where the expectation is the average over 7, see the last reference in 
[10]. This is a straightforward application of the usual rules of statistical 
mechanics 

This basic property, which we call the martingale property because 
in probabilistic jargon it would be translated as "(0)d\7]„;,] is a closed 
martingale for any observable O" is expected to survive in the continuum 
limit. In this limit, one has two powerful tools inherited from conformal 
invariance at hand, CFT to compute correlators and SLE to average 
over the piece of interface. Doing this for arbitrary values of k and c 
means mixing the degrees of freedom from two different models and there 
is a priori no reason for the martingale property to hold. An explicit 
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computation shows that it holds only if 

2kc= (6-k)(3k-8). 

The martingale property also gives information on the operator content 
of the CFT. As the discrete statistical mechanics examples show, there 
is a change in boundary conditions at the tip of the interface. The 
martingale property allows to identify this boundary changing condition 
operator as a primary field of weight h = {6 — k)/ {2k) degenerate at 
level two. 

This approach also exhibits a family of martingales which are at the 
heart of many probabilistic computations. For instance, it gives a sys- 
tematic way to interpret probabilities for SLE events as correlation func- 
tions of a CFT, and shows how the changes of behavior of the trace at 
K = 4, 8 are related to operator product expansions. The interested 
reader is refereed to [10]. 

2. LOEWNER CHAINS 

This section deals with more general 2D growth processes. Although, 
they do not fulfill the local growth and conformal invariance properties 
of SLEs, they are nevertheless described by dynamical conformal maps. 
We first present systems whose conformal maps have a time continuous 
evolution and give examples. We then go on by presenting a discrete 
version thereof in terms of iterated conformal maps. We explain inte- 
grability of Laplacian growth. The last part is a discussion concerning 
the limit of small ultraviolet cutoflF and the consequences of possible 
dendritic anomalies. 

2.1 CONTINUOUS LOEWNER CHAINS 

In this part, the exterior of the unit disk is used as the reference 
geometry. 

2.1.1 Radial Loewner chains. Let be a family of growing 
closed planar domains with the topology of a disk. Let = C \ be 

their complements in the complex plane. See Figure (1.16). To fix part 
of translation invariance we assume that the origin belongs to and 
the point at infinity to O^. 

Loewner chains describe the evolution of family of conformal maps 
ft uniformizing D = {w G C; \w\ > 1} onto O^. It thus describes the 
evolution of the physical domains Ot- We normalize the maps /j : D — >■ 
Ot by demanding that they fix the point at infinity, ft{oo) = oo and 
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that //(oo) > 0. With t parameterizing time, Loewner equation reads: 




The integration is over the unit circle {u £ C, \u\ = 1}. The Loewner 
density pt{u) codes for the time evolution. It may depends on the map 
ft in which case the growth process in non- linear. For the inverse maps 
gt = f^^ : Ot ^ T), Loewner equation reads: 

Q-^aAz) = -gt{z)f ^[-^^^^) Ptiu) (L8) 

The behavior of ft at infinity fixes a scale since at infinity, ftiw) ~ 
Rtw + 0(1) where Rt > 0, with the dimension of a [length] , is called 
the conformal radius of viewed from infinity. Rt may be used to 
analyze scaling behaviors, since Kobe 1/4-theorem (see e.g. [13]) ensures 
that Rt scales as the size of the domain. In particular, the (fractal) 
dimension D of the domains may be estimated by comparing their 
area At with their linear size measured by Rt'. At Rt' for large t - the 
proportionality factor contains a cutoff dependence which restores naive 
dimensional analysis. 




z - plane w - plane 

Figure 1.16 Uniformizing maps intertwining the 'physical' z-plane and the 'mathe- 
matical' ui-plane. 



2.1.2 The boundary curves. The boundary curve Tt = dOt is 
the image of the unit circle by ft- We may parameterize the boundary 
points by 'yt;a = ft{u) with u = e*". The Loewner equation codes for 
the evolution of the shape of and thus for the normal velocity of the 
boundary points. Only the normal velocity is relevant as the tangent 
velocity is parameterization dependent. The tangent to the curve is 
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r = iufl{u)/\fl.{u)\ and the outward normal is n = —it so that the 
normal velocity at 7^ is = ^e[ndtft{u)], or 

v,, = \n{u)\'!iic[dtftiu)/ufl{u)]. 

The r.h.s. is determined by the Loewner equation (1.7) because this 
equation may be viewed as providing the solution of a boundary value 
problem. Indeed, recall that for h{u) a real function on the unit circle, 
h{w) = § 2^ (^f^) h{u) is the unique function analytic outside the unit 

disk whose real part on the unit circle is h, i.e. 3fie[/t(?x)] = h{u). Thus, 
since dtft{w)/wfl{w) is analytic in D, the Loewner equation (1.7) is 
equivalent to: 

Vn = \fl{u)\pt{u) 

or more explicitly ^: 

{dalt;a) (dtTt^) - (dall^) {dt7t;a) = 2i |//(ii)|^ pt{u) (1.9) 

Hence, the evolution of the domain may be encoded either in the 
evolution law for its uniformizing conformal map as in eq.(1.7) or in 
the boundary normal velocity as in eq.(1.9). The two equations are 
equivalent. 

2.2 EXAMPLES 

2.2.1 Stochastic Loewner evolution (SLE). ^Prom this point 
of view, SLE is singular as it corresponds to a Dirac measure Sut (u) cen- 
tered at the position of a Brownian motion on the unit circle. The local- 
ity of this measure reflects the fact that SLE describes growing curves. 
Continuity of the Brownian motion reflects the absence of branching in 
the SLE growths. To make contact with the previous sections, the reader 
should check that eq.(1.5) becomes of the form (1.8) after the cylinder 
has been mapped to the outside of the unit disk. 

2.2.2 Laplacian growth (LG). This is a process in which the 
growth of the domain is governed by the solution of Laplace equation, i.e. 
by an harmonic function, in the exterior of the domain with appropriate 
boundary conditions. It originates from the hydrodynamical Hele-Shaw 
problem to be described below. 

To be precise, let P be the real solution of Laplace equation, V^P = 0, 
in Of with the boundary behavior P = — log + • • • at infinity and P = 



For SLE this equation has to be modified according to Ito calculus 



Loewner Chains. 29 



on the boundary curve Ft = dOt- The time evolution of the domain 
is then defined by demanding that the normal velocity of points on the 
boundary curve be equal to minus the gradient of P'. v^i — — (V-P)^^' 

This problem may be written as a Loewner chain since, as is well 
known, Laplace equation is solved via complex analysis by writing P as 
the real part of an analytic function. One first solves Laplace equation 
in the complement of the unit disk with the appropriate boundary con- 
ditions and then transports it back to the physical domain Ot using the 
map /(. This gives: 

P = -5fte $t with ^t{z) = log gt{z) 

The evolution equation for the map ft is derived using that the bound- 
ary normal velocity is f„ = — (VP)„. The above expression for P gives: 

Vn = -(VP)„ = \fi{u)\-' 

at point 7f = ft{u) on the boundary curve. As explained in the previous 
section, this is enough to determine dtft{w) for any \w\ > 1 since this 
data specifies the real part on the unit circle of the analytic function 
dtft{w)/wf[{w) on the complement of the unit disk. The result is: 

dtftiw) = Wf[{w) i fl,|2 (^) (1-10) 

J\u\ = l ^IT^U \jt{u)Y - 

It is a Loewner chain with pt{u) = \f[{u)\~'^. 

As we shall see below, Laplacian growth is an integrable system, which 
may be solved exactly, but it is ill-posed as the domain develops singu- 
larities (cusps ~ x^) in finite time. It thus needs to be regularized. 
There exist different ways of regularizing it. 

One may also formulate Laplacian growth using a language borrowed 
from electrostatics by imagining that the inner domain is a perfect con- 
ductor. Then F = 3f?e $t is the electric potential which vanishes on the 
conductor but with a charge at infinity. The electric field E = \/V 
\s E = Ex — iEy = dz^t- Its normal component E^ = |//('u)|~^ is 
proportional to the surface charge density. 

2.2.3 The Hele-Shaw problem (HS). This provides a hydro- 
dynamic regularization of Laplacian growth. The differences with Lapla- 
cian growth are in the boundary conditions which now involve a term 
proportional to the surface tension. It may be formulated as follows [14]. 

One imagines that the domain is filled with a non viscous fiuid, say 
air, and the domain Ot with a viscous one, say oil. Air is supposed to be 
injected at the origin and there is an oil drain at infinity. The pressure in 
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the air domain Kj is constant and set to zero by convention. In Ot the 
pressure satisfies the Laplace eqiiation V^P = with boundary behavior 
P = —(poo log |2;| + • • • at infinity reflecting the presence of the oil drain. 
The boundary conditions on the boundary curve are now P = —ant 
with a the surface tension and Kf the curvature of the boundary curve ^. 
The fluid velocity in the oil domain is u = —VP. Laplace equation 
for P is just a consequence of incompressibility. The evolution of the 
shape of the domain is specified by imposing that this relation holds on 
the boundary so that the boundary normal velocity is Vn = — (VP)n as 
in Laplacian growth. 

Compared to Laplacian growth, the only modification is the boundary 
condition on the boundary curve. This term prevents the formation of 
cusps with infinite curvature singularities. The parameter ^oo sets the 
scale of the velocity at infinity. In the following we set (p^o = 1. By 
dimensional analysis this implies that [time] scales as [length^] and the 
surface tension a has dimension of a [length]. It plays the role of an 
ultraviolet cut-off. 

A standard procedure [14] to solve the equations for the Hele-Shaw 
problem is by first determining the pressure using complex analysis and 
then computing the boundary normal velocity. By Laplace equation, 
the pressure is the real part of an analytic function, P = — Ke$t. The 
complex velocity v = Vx + iVy is v = dz^t- At infinity ^t{z) — log z + - ■ ■ 
and V 2± 1/z + ■ ■ ■. The boundary conditions on P demand that 



where i^tiw) is analytic in D with boundary value 3fte[t?t(?x)] = Kt{ft{u)) 
with Kt the curvature. Explicitly 



The evolution of ft is then found by evaluating the boundary normal 
velocity Vn = 3fie(V4')n at point 7t = ft{u): 



As above, this determines uniquely dtft{u) and it leads to a Loewner 
chain (1.7) with density: 



pt{u) = \fl{u)\-^(l + aet{u)) , et{u) = ^e[uduM'(^)] (1.11) 



^The curvature is defined by ft = —n.dsr/r'^ = Sm[rc3sT/|Tp] with r the tangent and n the 



{^t ° ft){w) = \ogw + (T-dtiw) 




= Mndz^t] = \ft{u)\-^ 5Re[l + auduOtiu)] 



normal vectors. An alternative formula is; t = ^SRe[l + 

R, the curvature is +1/R. 



]. For a disk of radius 
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The difference witli Laplacian growtfi is in tfie extra term proportional 
to a. It is highly non-linear and non-local. This problem is believed to 
be well defined at all times for a positive. 

2.2.4 Other regularized Laplacian growth (rLG). These reg- 

ularizations amount to introduce an UV cutoff 5 in the physical space 
by evaluating |//| at a finite distance away from dOt. A possible choice 
[20] is pt{uY/'^ = 5~Hrd{e : dist[/t(ii + ett); 30t] = 5). An estimation 
gives pt{u) X \f[{u + euu)\~'^ where goes to with 5, so that it naively 
approaches as (5 ^ 0. 

Another possible, but less physical, regularization consists in intro- 
ducing an UV cutoff u in the mathematical space so that pt{u) = 
\fl{u + vu)\-\ 

2.2.5 Dielectric breakdown and generalizations. A larger 
class of problems generalizing Laplacian growth have been introduced. 
Their Loewner measures are as in Laplacian growth but with a different 
exponent: 

pt{u) = \r^{u)\-^ , 0<a<2. 

Using an electrostatic interpretation of the harmonic potential, one usu- 
ally refers to the case a = 1 as a model of dielectric breakdown because 
the measure is then proportional to the local electric field En = |//(n)|~^. 
This is a phenomenological description. Just as the naive Laplacian 
growth these models are certainly ill-posed. They also require ultra- 
violet regularizations, one of which is described below using iterated 
conformal maps. 

2.3 CUSP SINGULARITIES IN LG 

The naive LG problem, without regularization, corresponds to the 
Loewner density pt{u) = \f[{u)\~'^. The occurrence of singularities may 
be grasped by looking for the evolution of domains with a Zn symmetry 
uniformized by the maps 

ftiw) = Rtw{i + ^^-n 

for some n > 2 and with < 1. This form of conformal maps is 
preserved by the dynamics. The conformal radius Rt and the coefficients 
f3t evolve with time according to dfRf = 2/(1 — /3f ) and Pt = {Rt/RcY~'^ 
with Rc some integration constant. The singularity appears when f3t 
touches the unit circle which arises at a finite time tc- At that time the 
conformal radius is Rc- 
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At tc the boundary curve T^^ has cusp singularities of the generic local 
form 



with ic a characteristic local length scale. In the present simple case 
ic ~ Re- At time t y tc, the dynamics is regular in the dimensionless 
parameter i~^\Jtc — t. The maximum curvature of the boundary curve 
scales as Kmax — ^c/ {tc — t) near tc and it is localized at a distance 
y^tT^^ away from the would be cusp tip. See Figure (1.17). 



This behavior is quite generic. As we shall see, conformal maps ft{w) 
such that their derivatives are polynomials in ■w~^ are stable by the 
Laplacian growth dynamics. By construction, their zeroes are localized 
inside the unit disk. A singularity in the boundary curve occurs if one 
of these zeroes converges to the unit circle. The singularity is then a 
cusp ic — as can be seen by expanding locally the conformal map 
around its singular point. 

Once regularized with an explicit ultraviolet cut-off, the processes are 
believed to be well defined for all time. The effect of the regularization 
procedure on the domain properties is presently unclear. The domain 
structures may a priori depend on how the problem has been regularized. 
In the hydrodynamic regularization -the Hele-Shaw problem- the cusp 
production is expected to be replaced by unlimited ramifications leading 
to dendritic growth. 

In the regularized model, the curvature of Fj is expected to remain 
finite at all time. Using scaling theory, a crude estimate of its maximum 
around the would be singularities may be obtained by interchanging the 
short distance scale \J {tc — t) near the singularity in the unregularized 
theory with the UV cutoff of the regularized theory. In the hydrody- 
namic regularization this gives Kmax — ^c/<7^ as a — > 0. 



Ic {5yf ~ {5xf 




Figure 1.17 Cups formation in Laplacian growth. 
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2.4 DISCRETE LOEWNER CHAINS 

2.4.1 DLA. DLA stands for diffusion limited aggregation [15]. It 
refers to processes in which the domains grow by aggregating diffusing 
particles. Namely, one imagines building up a domain by clustering 
particles one by one. These particles are released from the point at 
infinity, or uniformly from a large circle around infinity, and diffuse as 
random walkers. They will eventually hit the domain and the first time 
this happens they stick to it. By convention, time is incremented by 
unity each time a particle is added to the domain. Thus at each time 
step the area of the domain is increased by the physical size of the 
particle. The position at which the particle is added depends on the 
probability for a random walker to visit the boundary for the first time 
at this position. 

In a discrete approach one may imagine that the particles are tiny 
squares whose centers move on a square lattice whose edge lengths equal 
that of the particles, so that particles fill the lattice when they are glued 
together. The center of a particle moves as a random walker on the 
square lattice. The probability Q{x) that a particle visits a given site x of 
the lattice satisfies the lattice version of the Laplace equation V^Q = 0. 
It vanishes on the boundary of the domain, i.e. Q = on the boundary, 
because the probability for a particle to visit a point of the lattice already 
occupied, i.e. a point of the growing cluster, is zero. The local speed 
at which the domain is growing is proportional to the probability for a 
site next to the interface but on the outer domain to be visited. This 
probability is proportional to the discrete normal gradient of Q, since 
the visiting probability vanishes on the interface. So the local speed 
is Vn = (V(5)n- It is not so easy to make an unbiased simulation of 
DLA on the lattice. One of the reasons is that on the lattice there is 
no such simple boundary as a circle, for which the hitting distribution 
from infinity is uniform. The hitting distribution on the boundary of a 
square is not such a simple function. Another reason is that despite the 
fact that the symmetric random walk is recurrent is 2d, each walk takes 
many steps to glue to the growing domain. The typical time to generate 
a single sample of reasonable size with an acceptable bias is comparable 
to the time it takes to make enough statistics on loop-erased random 
walks or percolation to get the scaling exponent with two significant 
digits. Still this is a modest time, but it is enough to reveal the intricacy 
of the patterns that are formed. Fig. 1.18 is such a sample. The similarity 
with the sample in fig. 1.1, obtained by iteration of conformal maps, is 
striking. But a quantitative comparison of the two models is well out of 
analytic control and belongs to the realm of extensive simulations. 
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Figure 1.18 A DLA sample. 



DLA provides a discrete analogue of Laplacian growth. The particle 
size plays the role of an ultraviolet cutoff. Since DLA only assumes that 
the growth is governed by the diffusion of elementary patterns, DLA has 
been applied to a large variety of aggregation or deposition phenomena, 
see eg. [18]. 

During this process the clustering domain gets ramified and develops 
branches and fjords of various scales. The probability for a particle to 
stick on the cluster is much higher on the tip of the branches than deep 
inside the fjords. This property, relevant at all scales, is responsible for 
the fractal structure of the DLA clusters. 

Since its original presentation [15], DLA has been studied numerically 
quite extensively. There is now a consensus that the fractal dimension of 
2d DLA clusters is Ddia — 1-7L There is actually a debate on whether 
this dimension is geometry dependent but a recent study [19] seems to 
indicate that DLA clusters in a radial geometry and a channel geometry 
have identical fractal dimension. To add a new particle to the growing 
domain, a random walk has to wander around and the position at which 
it finally sticks is influenced by the whole domain. To rephrase this, for 
each new particle one has to solve the outer Laplace equation, a non- 
local problem, to know the sticking probability distribution. This is a 
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typical example when scale invariance is not expected to imply conformal 
invariance. 

2.4.2 Iterated conformal maps. As proposed in [17], an alter- 
native way to mimic the gluing of elementary particles consists in com- 
posing elementary conformal maps, each of which corresponds to adding 
an elementary particle to the domain. 

One starts with an elementary map corresponding to the gluing of a 
tiny bump, of linear size A, to the unit disk. A large variety of choices is 
possible, whose influence on the final structure of the domain is unclear. 
An example is given by the following formula {g\ is the inverse map of 
/a): 

, , z cos A — 1 

9x{z) = z - 

z — COS A 

fx{w)= = (2cosA)~^ W + 1 + Vw'^ - 2u;cos 2X + 1 

where f\ correspond to the deformation of the unit disk obtained by 
gluing a semi-disk centered at point 1 and whose two intersecting points 
with the unit circle define a cone of angle 2A. For A <C 1, the area of 
the added bump is of order A^. But other choices are possible and have 
been used. 

Gluing a bump around point e*^ on the unit circle is obtained by 
rotating these maps. The uniformizing maps are then 

fx-fiH = hi^^~'') 

The growth of the domain is obtained by successively iterating the 
maps /a„;6i„ with various values for the size A^ and the position 6n of the 
bumps. See Figure (1.19). Namely, if after n iterations the complement 
of the unit disk is uniformized into the complement of the domain by 
the map then at the next {n + 1)*'' iteration the uniformizing 

map is given by: 

Fin+1)H = i^(n)(/A„+,;(?„+iH) (1.12) 

For the inverse maps, this becomes 0(^^+1) — 9\n+r,0n+i ° ^(n)- 

To fully define the model one has to specify the choice of the parameter 
A„ and On at each iteration. Since A„ codes for the linear size of the 
added bump and since locally conformal maps act as dilatations, the 
usual choice is to rescale A„+i by a power of (e*^")| as: 

A„+i = Ao |i^„)(e^^")r"/2, 0<a<2 
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Figure 1.19 Iteration of conformal maps. 



The case q = 2 corresponds to DLA as the physical area of the added 
bump are approximatively constant and equal to Aq at each iterations. 
In the other case, the area of the added bump scales as 

The positions of the added bump are usually taken uniformly dis- 
tributed on the mathematical unit circle with a measure dO /2'k. 

It is clear that this discrete model with a = 2 provides a regularization 
of Laplacian growth with Aq playing the role of an ultraviolet cutoff. This 
may also be seen by looking at the naive limit of a small cutoff. Indeed, 
a naive expansion as A„ <C 1 gives that -F(n+i) = -F(n) + (^-^(n) with 

Cyj _|_ ^idn \ 
^^73^ ) 

where we used the expression of for A ^ 1. Using the recursive 
expression for A„ and averaging over 9 with a uniform distribution yields: 



{5F{w)) = XowF'iw) i —\F'{e 



27r w — e™ 

For a = 2 this reproduces the Loewner chain for Laplacian growth. But 
this computation is too naive as the small cutoff limit is not smooth, a 
fact which is at the origin of the non trivial fractal dimensions of the 
growing domains. 

There are only very few mathematical results on these discrete mod- 
els. The most recent one [27] deals with the simplest (yet interesting 
but not very physical) model with a = 0. It proves the convergence 
of the iteration to well-defined random maps uniformizing domains of 
Hausdorff dimension 1. However, these models have been studied nu- 
merically extensively. There exists a huge literature on this subject but 
see ref.[24] for instance. These studies confirm that the fractal dimension 
of DLA clusters with a = 2 is -Ddia — 1-71 but they also provide further 
informations on the harmonic measure multi-fractal spectrum. Results 
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on the a dependence of the fractal dimension may be found eg. in ref. 

Various generahzations have been introduced. For instance, in ref. [26] 
a model of iterated conformal maps has been defined in which particles 
are not added one by one but by layers. These models have one control 
parameter coding for the degree of coverage of the layer at each iterative 
step. By varying this parameter the model interpolates between discrete 
DLA and a discrete version of the Hele-Shaw problem. The fractal di- 
mension of the resulting clusters varies with this parameter and seems 
to indicate that the fractal dimension of this discrete analogue of the 
Hele-Shaw problem is 2, a point that we shall discuss further in the last 
section. 

2.5 INTEGRABILITY OF LAPLACIAN 
GROWTH 

Laplacian growth is an integrable system, at least up to the cusp 
formation. Let us recall that it corresponds to a Loewner chain with a 
density pt{u) = \ft{u)\~^, or equivalently to the quadratic equation 

{dalt;a) {dtlt^) - (dajt^) {dtjt;a) = 2^ (1.13) 

for the dynamics of the boundary points ^yt-a = fti^), u = e*". What 
makes the model integrable is the fact that the r.h.s of eq.(1.13) is con- 
stant. Eq.(1.13) is then similar to a quadratic Hirota equation. Hints on 
the integrable structure were found in [21] and much further developed 
in [22, 23]. 

2.5.1 Conserved quantities. We now define an infinite set of 
quantities which are conserved in the naive unregularized LG problem. 
They reflect its integr ability. We follow ref. [22, 23]. These quantities 
may be defined via a Riemann-Hilbert problem on Ft specified by, 

5+(7)--S_(7) = 7 , 7ert (1.14) 

for functions S- and S-^. respectively analytic in the outer domain Ot and 
in the inner domain Kj. We fix normalization by demanding S- (oo) = 0. 
We assume Tf regular enough for this Ricmann-Hilbcrt problem to be 
well defined. As usual, S± may be presented as contour integrals: 



Tt 2z7r z — 7 



The conserved quantities are going to be expressed in terms of S±. We 
thus need their time evolution. Differentiation of eq.(1.14) with respect 
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to time and use of the evolution equation (1.13) gives: 

dtS+{-f)-dtS.{j)=2g'^{j)/gt{j) 

Notice now that g'til) / dtil) is the boundary value of (loggt)' which by 
construction is analytic in Of. We may thus rewrite this equation as a 
trivial Ricmann-Hilbert problem, dtS+{'j) — {dfS- + 2(Iogg't)')(7) = 0) 
so that both terms vanish: 

dtS+{z) = and {dtS- + 2{loggty){z) = (1.15) 

Since 5*+ is analytic around the origin, we may expand it in power 
of z. Equation ^tS-^-{z) = then tells us that S+{z) is a generating 
function of conserved quantities: S+{z) = J2k>o ^^-^k with 

4 = £ 1^77"'"' ' 9th = 0. (1.16) 

This provides an infinite set of conserved quantities. 

Since S- is analytic around infinity, it may be expanded in power of 

l/z: S-{z) = —At/TTZ H with At = ^^77 °f ^t- '^^^ 

second equation {dtS-+2{\og gty){z) = with gt{z) = R^^z+0{1) then 
implies dtAt = 27r. The area of the domain grows linearly with time, up 
to the time at which the first cusp singularity appears. 

2.5.2 Simple solutions. A particularly simple class of conformal 
maps, solutions of the Laplacian growth equation, are those such that 
their derivatives are polynomials in w~^. They may be expanded as: 

N 

ft{w) = Y.fnW^'"^ fo = Rt>0 (1.17) 

n=0 

with A'^ finite but arbitrary. The dynamical variables are the A + 1 
coefficient fo, - ■ ■ ,fN- They are all complex except /o which is real. It 
will be convenient to define the function ft by ft{w) = ft(w). 

The fact that this class is stable under the dynamics follows from the 
Loewner equation (1.10). The trick consists in using the fact that the 
integration contour is on the unit circle so that \ft{u)\^ = fl{u)fl{l/u). 
The contour integral then involves a meromorphic function of u so that 
it can be evaluated by deforming the contour to pick the residues. This 
is enough to prove that dtft{w) possesses the same structure as ft{w) 
itself so that the class of functions (1.17) is stable under the dynamics. 

Alternatively one may expand the quadratic equation (1.13) to get a 
hierarchy of equations: 

XI (1 - ")[/n/j+n + fnf-j+n] = 2(5j;0 
n>0 
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For j = 0, this equation tells us again that the area of the domain grows 
linearly with time. Besides this relation there are only N independent 
complex equations for j = 1, • • • , AT which actually code for the conserved 
quantities. 

To really have an integrable system we need to have as many inde- 
pendent integrals of motion as dynamical variables. Thus we need to 
have conserved quantities. These are given by the I^'s defined above 
which may rewritten as 

r du f't{u)ft{i/u) 

^ J\u\=i2m ft{uf+^ 

Only the first N quantities, Iq, - ■ ■ ,In-i = R^''^ /n are non-vanishing. 
They are independent. They can be used to express algebraically all 
/n's, n > 0, in terms of the real parameter /o = Rt- The area law, 

At = 7r[i?f + 5^(1 - n)\fnf] = 2Trt, 

n>l 

with the fnS expressed in terms of Rt, then reintroduces the time vari- 
able by giving its relation with the conformal radius. 

2.5.3 Algebraic curves. As was pointed out in [23], solutions of 
Laplacian growth and their cusp formations have an elegant geometrical 
interpretation involving Riemann surfaces. 

Recall that given a sufficiently smooth real curve Tt drawn on the com- 
plex plane one may define a function S{z), called the Schwarz function, 
analytic in a ribbon enveloping the curve such that 

'5(7) =7, 7ert 

By construction, the Schwarz function may be expressed in terms of 
uniformizing maps of the domain bounded by the curves as S{z) = 

UVgtiz)). 

The Riemann- Hilbert problem (1.14) defining the conserved charges 
then possesses a very simple interpretation: S± are the polar part of the 
Schwarz function respectively analytic inside or outside Ff, i.e. S{z) = 
Sj^{z) — S-{z). Thus the polar part 5*+, analytic in the inner domain, 
is conserved. The polar part S'_, analytic in the outer domain, evolves 
according to eqs.(1.15). Since \oggt{z) is analytic in the outer domain, 
these equations are equivalent to the single equation: 

dtS{z) = -2{\oggt{z))' (1.18) 

Now the physical curve F^ may be viewed as a real slice of a complex 
curve, alias a Riemann surface. The latter is constructed using the 
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Schwarz function as follows. Recall that s = S{z) is implicitly defined 
by the relations z = ft{w), s = ft{l/w). In the case of polynomial 
uniformizing maps we get the pair of equations 

Z = foW + fl + f2W-^ + ■ ■ ■ + /nW^-"" 
S = f0W-^ + fl + f2W + ---+fNw''-^ 

Eliminating w yields an algebraic equation for z and s only: 

R: R{z,s) = (1.19) 

with R a polynomial of degree N in both variables, R{z, s) = Jnz^ + 
Ins^ + • • •• Eq.(1.19) defines an algebraic curve R. It is of genus zero 
since by construction it is uniformized by points w of the complex sphere. 
It has many singularities which have to be resolved to recover a smooth 
complex manifold. 

The Riemann surface R may be viewed as a A^^-sheeted covering of the 
complex z plane: each sheet corresponds to a determination of s above 
point z. At infinity, the physical sheet corresponds to ^ ~ Jqw with 

— > oo so that s ~ {z/ fo)^^^ /at, the other N — 1 sheets are ramified 
and correspond to 2; ~ In/w^^^ and s ~ fo/w with w — ^ so that 
z ~ {s/ fo)^~^ Jn- Hence infinity is a branch point of order — 1. 

By the Riemann-Hurwitz formula the genus g is 2g — 2 = —2N + u 
with 1/ the branching index of the covering. Since the point at infinity 
counts for i^oo = N—2, there should be N other branch points generically 
of order two. By definition they are determined by solving the equations 
R{z, s) = and dsR{z, s) = 0. Since the curve is uniformized hy w & C, 
these two equations imply that z'{w)dzR{z{w), s{w)) = 0. Hence either 
z'lw) = 0, dzR 7^ 0, and the point is a branch point, or z'{w) 7^ 0, 
dzR = = dsR, and the point is actually a singular point which needs 
to be desingularized. So the N branch points at finite distance are the 
critical points of the uniformizing map z = ft{w). 

The curve R possesses an involution (z, s) {s, z) since R{s, z) = 
R{z, s) by construction. The set of points fixed by this involution has 
two components: (i) a continuous one parametrized by points w = u, 
\u\ = 1 -this is the real curve Ft that we started with- and (ii) a set of 
A'^ isolated points which are actually singular points. 

The cusp singularity of the real curve arises when a isolated real 
point merges with the continuous real slice Tf. Locally the behavior is 
as for the curve = ev'^ + v"^ with e ^ 0. 

The simplest example is for N = 3 with Z3 symmetry so that ft{w) = 
w + b/w^ and 

w'^ z = + b , w s = 1 + bw^ 
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We set /o = 1 and /s = b. Without lost of generality we assume b real. 
The algebraic curve is then 

R{z, s) = bz^ + bs^ - b^s^z^ + (6^ - i){2b^ + i)sz - {b^ - 1)^ = 

Infinity is a branch point of order two. The three other branch points 
are at z = 3lo (6/4)^/^, s = lo^ (26^ + 1)(26)~^/^ corresponding to w; = 
(26)1/3 with LO a third root of unity. They are critical points of z{w). 
There are three singular points at z = (1 — b'^)/b, s = a;^ (1 — b'^)/b 
corresponding to it; = uj{l it \/l — 46^)/26. The physical regime is for 
b < 1/2 in which case the real slice Tt = {z{u)., \u\ = 1} is a simple 
curve. The singular points are then in the outer domain and the branch 
points in the inner domain. The cusp singularities arise for b = 1/2. For 
6 > 1/2 there are no isolated singular points, they are all localized on 
the real slice so that Vt possesses double points. See Figure (1.20). 




Figure 1.20 Subcritical and critical algebraic curves. Black circles are singular 
points. Crosses are branch points. 



2.6 DENDRITIC ANOMALIES 

We now would like to discuss a few points concerning the regularized 
models and their small cut-off limits -which arc actually what we are 
interested in. We shall point out that the hydrodynamic regularization 
used in the Hele-Shaw problem possesses essential differences with that 
used in other regularized models, say DLA. This opens the possibility 
for the Hele-Shaw problem and DLA not to be in the same universality 
class. 

This observation is based on the conjectural existence of anomalies 
in the Hele-Shaw problem -a word which refers to quantities which, 
although naively vanishing in the small cut-off limit, are actually non 
zero but finite in this limit due to compensating effects. 
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An analogy with Burgers turbulence may be useful. The Id Burgers 
equation is dtu + udxU — ud^u = for a velocity field u{t,x). The 
viscosity v plays the role of ultraviolet cut-off. At = this equation is 
simply Euler equation which may be easily solved. At = 0, any smooth 
initial data produces shocks in finite time with a discontinuity of u. In 
the regularized equation with i/ 7^ 0, there is no shock production but 
the field e = u{dxu)'^^ which naively vanishes when 1/ ^ 0, has actually 
a non trivial limit. It does not vanish because there is a compensation 
between small viscosity v and large velocity gradient (dxu). The field e 
codes for the amount of energy dissipated inside the shocks. 

The analogy that we would like to put forward is that the cusp for- 
mation in Laplacian growth is analogous to shock formation in Burgers 
and that the surface tension a in Laplacian growth plays the role of the 
viscosity in Burgers. There are potential compensating effects between 
small surface tension and large curvature in the Hele-Shaw problem. 

2.6.1 Almost conserved quantities. Wc now describe the ef- 
fect of the regularization procedure on the conserved quantities of Lapla- 
cian growth. Recall that they were defined starting from the Riemann- 
Hilbert problem: 

/i+(7)-/i-(7) = 7 , TeTi (1.20) 

for functions /i_ and respectively analytic in and K^. Wc use the 
notation h± instead of S± in order to avoid confusion with the unreg- 
ularized Laplacian growth. We again fix normalization by demanding 
h-{oo) = so that h±{z) = — /p^ ^ Time differentiating eq.(1.20) 
using the evolution equation (1.9) now gives: 

with Vn the boundary normal velocity -in LG, we had Vn = If/til)]- 

In the case of the Hele-Shaw (HS) problem, the hydrodynamic regu- 
larization of Laplacian growth, the normal velocity possesses a particular 
form (1.11) involving the derivative of 1?. Eq.(1.21) may then be rewrit- 
ten in terms of the potential, ^t{z) = log gt{z) + a-dt{gt{z)), and of the 
curvature as: 

dth+{j) - dth-ij) = 2$;(7) - 2aK[{j) 



Since is analytic in Ot, the solution of this equation is given by that of 
the Riemann-Hilbert problem for Kt- So, let k- and k+ be respectively 
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analytic in Ot and Kt with boundary values k+ — /c_ = Kt on Tf. Then: 



These equations codes for the behavior of h± in the limit of vanishing 
surface tension. 

A noticeable point is that the above r.h.s. are total derivatives. In 
particular, contribution to the 1/z term in dth- only comes from the 
logarithmic factor in the potential. This implies that at a ^ the area 
of the domains also grows linearly with time, At — ^irt at large time ^, 
and it is independent of the UV cutoff a. This is contrast with other 
regularizations for which the area growth may be cutoff dependent. 

2.6.2 A conjecture. The supplementary terms proportional to 
(7Kj_ in the quasi conservation laws (1.22) naively vanish when the UV 

cutoff is removed. Dendritic anomalies refer to the possibilities for these 
terms not to vanish as the cutoff goes to zero, say a — > in the Hele- 
Shaw (HS) problem. Such possibility arises due to a balance between 
small surface tension and large curvature. 

Recall the previous HS estimation of the maximum curvature /^max — 
^c/c^ with ic the local scale associated to the cusp. These maxima are 
localized around points 7c which would turn into cusp singularities in 
absence of surface tension. So, we may expect that in the HS problem 
the product <tk(7) is a function picked around its maximum values and 
is of the form ^(p{^^^^^) when the point 7 moves away from jc with ip 
finite at the origin and decreasing rapidly away from it. For isolated 
singular points 7c, this would converge in the limit cr to a sum of 
Dirac point measures. 

This leads us to conjecture that the anomalous term (Tk(7) goes to 
a non trivial but finite distribution as o" ^ 0. This distribution then 
contributes to eq.(1.22) to break the 'classical' conservation laws valid 
in the unregularized theory, hence the name dendritic anomaly. 

The Loewner equation (1.7) -or the evolution equation (1.9)- are 
then expected to have a finite limit as o" ^ with a non trivial right 
hand side. The solution would be well-defined for all times and would 
describe a non trivial familly such that, for fixed t, is the limit of 
the domains as cr ^ 0. This is consistent with the observation that 
the area of is independant of a and grows linearly with time for any 



^This result concerning the area is simply a consequence of the fluid incompressibility. But 
the general method allows to deal with the complete hierarchy of quasi conserved charges. 



dth+{z) 
dth.{z) 



—2adz 
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£7, ie. At — 2Trt, so that the existence of a hmiting set is not ruled 
out. 

It is implicit in the previous conjecture that the uniformizing maps 
ft have a finite limit, say ft, as cr ^ such that the image of the unit 
circle by ft is the boundary of C \ K^. This leads us to conjecture that 
the limiting domain Of, which is the component of C \ containing 
infinity, is uniformized onto D by = 

Let E4 be the conformal radius computed at finite a using the maps ft- 
Recall that if Dhs is the fractal dimension we have i?^''^ x At for large 
t. By dimensional analysis Rt alAt/cr"^)^^^^^ since a has dimension 
of a [length] . Because in the HS problem the area grows linearly with 
time without any dependence on a, this also reads: 

Since we expect the limiting domains and the limiting conformal maps 
to exit, Rt should have a finite limit as a — > which coincides with the 
conformal radius Rt computed using ft. As consequence, we expect: 

^hs = 2 

That is: the fractal dimension of Laplacian growth clusters, within the 
hydrodynamic rcgularization, is 2. If these conjectures are true, they 
imply that the Hele-Shaw problem is not in the same universality class 
as DLA. 

These conjectures are compatible with that of ref.[26] based on nu- 
merical studies of generalized iterated conformal maps. They are not 
expected to apply to arbitrary regularizations because in these cases the 
area may depend on the UV cutoff and the previous dimensional analy- 
sis does not apply. It would then be natural to wonder whether we may 
define renormalized uniformizing maps in the limit of a vanishing cutoff. 
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